Partial regularity for the Navier-Stokes-Fourier 

system 



L. Consiglieri* 
October 24, 2009 

Abstract 

This paper addresses a nonstationary flow of heat-conductive in- 
compressible Newtonian fluid with temperature-dependent viscosity 
coupled with linear heat transfer with advection and a viscous heat 
source term, under homogeneous Dirichlet boundary conditions. The 
partial regularity for the velocity of the fluid is proved to each proper 
weak solution, that is, for such weak solutions which satisfy some lo- 
cal energy estimates in a similar way to the suitable weak solutions of 
the Navier-Stokes system. Finally, we study the nature of the set of 
points in space and time upon which proper weak solutions could be 
singular. 
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1 Introduction 

Due to diverse applications, coupled systems governing incompressible flows 
are subject of intensive analytical and numerical investigation (see for in- 
stance El El dH [El HH1 EQl EH EH] and the references therein). Here we 
deal with a nonstationary flow of heat-conductive incompressible Newtonian 
fluid with temperature-dependent viscosity coupled with linear heat transfer 
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with advection and a viscous heat source term, under homogeneous Dirichlet 
boundary conditions, and we study the partial regularity for proper weak 
solutions to the coupled system under study. The proper weak solution is 
each weak solution that satisfies some local energy estimates in a similar way 
to the suitable weak solutions of the Navier-Stokes system [8] . Although the 
techniques used in the present work could be considered standard, the result 
is new because the mentioned techniques can not be directly applied. Un- 
der isothermal effects, the viscosity is constant and the problem is described 
by the Navier-Stokes equations. We refer to [22j [23J, [21] where the par- 
tial regularity theory is studied: at the first instant of time when a viscous 
incompressible fluid flow with finite kinetic energy in three space becomes 
singular, the singularities in space are concentrated on a closed set whose 
one dimensional Hausdorff measure is finite. This investigation character- 
izes some geometric properties and measures theoretic properties of the sets 
of points in space and time upon which weak solutions could be singular 
[3], [16]. In such conditions, it is known that the concept of weak solutions in 
the sense given by Leray and Hopf is not sufficient to establish their partial 
regularity. Different regularity criteria for suitable weak solutions to the N-S 
system has been introduced in terms of the smallness of functionals that are 
invariant with respect to the natural scaling either the velocity or its gradient 
[TUl \13\ [T8] . We prove the existence of regular points in the sense due to [15] , 
that is, the function is Holder continuous on a parabolic cylinder centred at 
such point. We remark that in the popular definition given at [5] the Holder 
space is replaced by the space of essentially bounded functions. 
We recall that the following continuous inclusion 

K\Q T ) ■= {v E L«(0,T;W*«(fi)) : d t v E L q (Q T )} - C k ^{Q T ) 

only occurs if g > (n+2)/(2 — k). This means for k = that q > n/2 + 1 > 2 
(n = 2,3), i.e., the Banach space (Qt) is not embedded in the Banach 
space of Holder continuous functions with exponent a in the x-variables and 
a/2 in the t- variable. It is known that the boundedness of the velocity of 
the fluid, when a weak solution of N-S equations is care of, requires the 
so-called Ladyzhenskaya-Prodi-Serrin sufficient condition of the norm of the 
mixed Lebesgue space L s,t {Qt) with 2/r + n/s = 1 if s > n. Likewise 
recently a supplementary criterion concerning the N-S equations ensures the 
boundedness of solutions [2], representing the pressure by a relation on the 
velocity which envolves its spatial derivatives of second order. We emphasize 
that even in 2D, the regularity theory for the second derivative in space 
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or the derivative with respect to the time variable is not available to N-S- 
F system because the viscosity is a non-H61derian function. It remains to 
know whether or not it is possible to find a L°°-bound of the velocity of 
heat-conducting fluids. 

The partial regularity up to the boundary for the N-S-F system is still 
another open problem, since the study for N-S equations is provided by a 
constant viscosity in order to request that the spatial derivative of second 
order, V 2 w, belongs to a convenient mixed Lebesgue space, namely Ls>2(Q T ) 
[26l [27] what in the present study does not happen. Even in the study of 
parabolic equations the partial regularity up to the boundary is proved by 
means of the existence of d t u G L 2 {Q T ) PQ. Once more such regularity does 
not occur in the present study. 

Let Q C W 1 be a bounded open domain, dQ its boundary, and T > 0. 
Let us consider the initial boundary value problem of the N-S-F system: 

d t u - div(/i(0)Z>u) + (u - V)u = f - Vp in Q T := fix]0, T[; 

div u = in Q T ; ^ ' 

d t 9 - kA9 + u • V0 = fi{9) \Du\ 2 in Q T ; (2) 

u lt=o = u o> 0\ t =o = °o in fi; (3) 
u = 0, 9 = ondQx]0,T[, (4) 

where 9 represents the temperature, u the velocity of the fluid and Du = 
|(Vu + Vu T ), p denotes the pressure, \i the viscosity, f denotes the given 
external body forces, k denotes the conductivity assumed to be a fixed pos- 
itive constant, and u and 9 Q are some given functions. For simplicity, the 
density and the heat capacity are constants assumed equal to one and we 
do not consider the existence of the external heat source, since the heating 
dissipative term is the main mathematical difficulty The product of two 
tensors is given by D : r = DijTij (in Einstein's convention) and the norm 
by \D\ 2 = D : D. 

The outline of the paper is as follows. Next section is concerned to the 
presentation of the appropriate functional framework and the main results. 
The Section H] is devoted to the proof of some auxiliary results introduced at 
Section [3j In Section [51 we prove the partial regularity result (Theorem 12. 2j) . 
Finally an estimate to the parabolic Hausdorff dimension on the set of the 
singularities for the fluid velocity is provided in Section [HI 
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2 Assumptions and main results 



Here we assume that Q C M n (n = 2, 3) is a bounded open domain sufficiently 
regular, e.g. of class C 2 . For 1 < q < oo, we introduce the Sobolev space of 
divergence-free fields 

Jl' q (tt) = {u G Wj'^fi) : V ■ u = in fi} 

with norm 

II ' ||l,g,fi = II V • \\q t Cl- 

We denote by bold the vector spaces of vector-valued or tensor-valued func- 
tions. For any set A, we write (u,v)a_ '■= Ja uv whenever u G L q (A) and 
v G L q '(A), where q' = q/(q — 1) is the conjugate exponent to q, or simply 
(-, •) whenever there exists no confusing at all, and we use the symbol (-, •) to 
denote a generic duality pairing, not distinguished between scalar and vector 
fields. 

Definition 2.1 We say that the triple (u,p,8) is a weak solution to the 
Navier- Stokes- Fourier (N-S-F) problem (TjP-(gp in Q T if p G L^ n+2 ^ n (Q T ), 

uGW:= L°°(0, T; L 2 (ft)) n L 2 (0, T; Jo' 2 (fi)), 
6 G S := L°°(0, T; L\Q)) n L«(0, T; g < 

t u EX := L 2 (0,T; W -1 ' 2 (fi)) n L (n+2)/n (0, T; W' Un+2)/n (Q)), 



£ q(n + l) 2(n + 2)' 
and satisfies the variational formulation 

(d t u, v) + / (^(fl)ZJu :Dv+(u-V)u-v) dxdt = 

= / (f • v + p div \)dxdt (5) 

Vv G L°°(0, T; Wj ,0O (fi)), u| t=0 = u in fi; 
(d t 6,(f))+ I (kV6-9u) -V(pdxdt= [ fi(9)\Du\ 2 (f)dxdt (6) 

V0 G L°°(0, T; W 1>QO (fi)), 0| t=o = # in fi. 
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In (J5]) the convective term verifies (u • V)u G 'h( n+2 ^ t < n+1 \Q T ), for every 
u G U ^ L 2 ( n+2 )/ n (g T ). In ® the advection term 6u G L e (Q T ) if 9 G 
g ^ L<?("+i)/«(Q T ) for £ > 1, i.e. q > 1 if n = 2 and q > 15/14 if n = 3 
(q > 2n(n + 2) / ((n + 4)(n + 1))). For n = 2, 3, we remark that the embedding 
holds X ^ L( n+2 )/ n (0,T; W- 1 -( n+2 )/ n (fi)). In conclusion, all terms in the 
above equalities are meaningful. 

We begin by recalling the existence of weak solutions such that verify 
some local energy inequalities (the proof of the existence result can be found 
in [8]). We assume that the following hypotheses hold 

(Al) f : Q T -y W 1 is given such that f G L 2 ( 1+e °)(Q T ) with e > 0; 

(A2) fi : K. — > R is a continuous function such that 

< fi# < n(s) < fi*, Vs G E; (7) 

(A3) u G L 2 (fi), 6 G L x (p) such that 

div uo = in Q; ess inf 9q(x) > 0. (8) 



Theorem 2.1 Under the assumptions (A1)-(A3), the N-S-F problem defined 
by (EP- has proper weak solutions in the following sense. The weak solution 
in accordance to Definition \2. 1\ satisfies the following local energy inequalities 



I |u(:M) a l 2 ^ (rM)&+ I ^(6)\Du\ 2 V 2 dxdT< 
Jn 2 JQ t =Qx]o,t[ 

<2 fi(9)ipDu : ((u-a) (g) V(p) + / |u - a|V($P + u ■ V<£>) + 
JQt JQt 

+2 / ptpu ■ V^pdxdr + f • (u — &)(p 2 dxdT, (9) 
JQt JQt 



(C + # 2 ) 3 / 2 



< / + # 2 <M + kA4> + u • W) + / /i(0)|£u 



2 ^ (io) 



(C + # 2 ) 1/2 



/or a// G C$°(Q T ), for any a G R n , /or any ( > 0, and for all ip G C™(Q T ) 
such that ip > 0, a.et g]0,T[. Moreover, 6 > in Qt- 
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REMARK 2.1 The local energy inequality for the temperature, / T70j) . plays an 
essential role in the proof of the decay lemma (cf. Theorem VJ. 1\) because if 
we take C -> 0+ we obtain 6 G L°°(0, T; L 1 1 oc (fi)). 

In order to establish the partial regularity, let us introduce some addi- 
tional notations. 

Definition 2.2 Given a positive number X, we set 



c\{f,u) := sup{ 




R>0, Q(z,R) CCtu} (11) 



for all u C Qt- In particular, we set c\(f) = c\(f, Qt)- 

Given a positive number A we introduce the "parabolic" variant of the 
Morrey spaces [15] 

M 2 , X (Q T ) = {f G Ll c (Q T ) : c x (f,Q T ) < oo}. 

Let Zq = (x , t ) G Q T and R > be such that Q(z , R) = B(x , R)x]t — 
R 2 ,to[c Qt- For any set A, we denote by \A\ the usual Lebesgue measure 
of the set, and we write -f A v := Aj Ja v wnenever v £ ^(A), or simply 
{v)a — -f* v - Observing that 2(n + 2)/(n + 4) < a' < g(n + l)/n is equivalent 
to 

g(w + 1) 2(n + 2) _ f 4 if n = 2 

g(n+l)-n n \ 10/3 if n = 3 

we can choose 1 < g < (n + 2)/(n + 1) such that q > 3n/(2(n + 1)), i.e. 
g(n + l)/(g(n + 1) — n) < 3. So we set 

3 < a < 2(ra + 2)/ra. (12) 

Let us introduce the following scaling invariant functional: 

Y R (z ;u,p,9) = (j QM \u\ a dz} ' +R( y jQ {z(hR) \p\ n ¥dzy +2 + 

The choice of the above exponents is consequence of u 6 h 2( - n+2 ^ n (QT), 
p G L( n+2 )/ n (Q T ), 6 G L«( n+1 )/ n (Q T ) and ufl G L x (g T ). We refer [Tg] in order 
to compare to the those exponents in the analysis for the N-S system into 
the three-dimensional space (n = 3). 
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Theorem 2.2 Suppose A and a are positive numbers such that |Hj) is ver- 
ified. Let f G M 2j \(Qt) and (Al)-(AS) be fulfilled. Assume that (u,p,9) 
is a proper weak solution of the N-S-F system in Qt- Then, there exists a 
positive constant A depending only on , A andn. Moreover, if for 

any z G Qt 

lim sup RY R (z ; u,p, 9) < A, (13) 

R-+0+ 

then z Q is a regular point in the following sense: the function z i— > u(z) is 
Holder continuous in some neighbourhood of the point zq. 

The local energy inequality for the temperature has a nonstandard format. 
It is still an open problem the proof that the pointwise criterion 

limsup^- ( [ |Vu| 2 rf^ +( f \Ve\ q dz] < 7* 

r^O R \J Q ( ZOtR ) J \JQ(z ,R) J 

implies the local condition for the regularity ([TBI , for some constant 7* > 0. 

Remark 2.2 Under all conditions of Theorem \2.2\ we denote by Qo the set 
of all regular points such that satisfy By definition, Qo is an open set. 

The set S = Qt \ Qo, known as the singular set, reads 

S = {z G Qt '■ lim sup RYr(z; u,p, 6) > A}. 

K->0+ 

Moreover, it can be proved that it is of Lebesgue measure zero. 

The parabolic Hausdorff dimension on the singular set can be estimated. 
Indeed, we state the following result. 

Theorem 2.3 Let a be given as in (TJ|j. Then 7fM<*- 2 )/(a-i) (5-) = 0, where 
ri d is the d- dimensional Hausdorff measure with respect to the standard 
parabolic metric of a set S defined as follows 

00 00 
H d (S) := liHiinf^.Rf : S C (jQ{zi,Ri), 0<Ri<5}. 

^° i=l i=l 

Therefore, the parabolic Hausdorff dimension of S does not exceed n — (a — 
2)/(a - 1), i.e. dim H (S) := inf{d > : H d (S) = 0} < n - (a - 2)/(a - 1). 
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3 Auxiliary results 

The main tool for the partial regularity analysis is the decay property of 
the scaled Lebesgue norms of the triple velocity-pressure-temperature which 
is based on a standard "blow up" method and the decomposition of the 
pressure. Let us introduce the following scaling invariant functional: 

/ \ l/a 

Y R (z ;u, Pl 6) = i Q{z0jR) \u(x,t) - (u) Q{zihR) \ a dzj + 

/ \ n/(n+2) 

+ R { J Q ( Z0 ,R) Ip(<M) - (p)B( X0 ,R)( t )\ (n+2)/ndz ) + 

+R{i Q{Z0 , R) \^t)-(e) Q(Z0 , R) fdzy /a . 

Theorem 3.1 (Decay estimate) Suppose that < q < 1 and < (3 < A 

are fixed numbers. For each function \x satisfying there exist positive 
constants 7 and e > depending only on fi#, yfi 1 , q, (3 and A such that for 
any proper weak solution (u,p,9) of the N-S-F problem in Qt verifying 

VQ(z,R) CC Q T : < R < e, ] 
R\(u) Q(ZtR) \ < 1, R\{9) Q{z , R) \< 1 /4 \ (14) 
Y R (z;u,p,9) + c x (f)RP <^ J 

we have the decay estimate 

Y iR (z\ u,p, 6) < C*q a (Y R (z; u,p, 6) + c x (f)R p ) 

with < a < 2 — n/2 and some absolute constant C* depending only on //# 
and //#. 

To prove Theorem 12.21 some iterative estimates are required [T5l Lemmas 
2.5, 3.1]. Let us state an iterative result, in which a careful choice of q 
significantly simplifies the iterative formula as well as its proof (compare to 
[IS Lemmas 2.5, 3.1]). 

Lemma 3.1 Under all conditions of Theorem \3.1\ let < 5 < a. If addi- 
tionally q is such that 

C^ a - 5 < 1/2 and + 1, (15) 

then for any proper weak solution (u,p, 6) of the N-S-F system in Qt satis- 
fying (73) verify 
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1. for any m G N, 

Y, mR (z;u,p,6) <<; m5 (Y R (z;u,p,6) + c x (f)R 13 ); (16) 

2. for allpe]Q,qR], 

Y p (z; u,p, 6)<C (£) (Y R (z; u,p, 6) + c A (f)i^), (17) 

where C = C(s,5) denotes a positive constant. 

Finally, we state the additional technical result. 

Lemma 3.2 Assume that (u,p,9) satisfies the system in the sense 

of distributions. For every R, r > 0, e = (yo,s ) and z = (x ,to) suc h that 
Q(z ,R) C Q T , 

1. if we introduce the functions 

u R {y,s) = ^ u {x + f{y-y ),t +(f) 2 {s-s )), 
p R (y, s) = (f ) 2 p(x + f (y - y ), *o + (f)' ( s ~ s o)), 
9 R (y, s) ee (f ) 2 ^(x + f (y - y ), to + (f ) 2 (a - *„)), 

i/ien (u R ,p R ,6 R ) satisfy the transported system in Q(eo,r) 

d s u R - div s/ (/i R (^) J D y u R ) + (u R • V y )u R = f fl - V yP R ; (18) 

div y ii fl = 0; (19) 

d s 9 R - kA y 6 R + u R ■ V y 6 R = /2 R (e R )\D y u R \ 2 ; (20) 

with the function /^(i?) ee fi{r 2, &/ R 2 ) and 

f R (y,s)=(^\ f( Xo + ^-y ) )to +(j) (s-so)); 

ry r (e ;u fl ,p ,l ,O = i?? fi (z ; u,p, #). (21) 
Moreover ^ ^ 

c x (f R ,Q(e ,r))< (^j c x (f,Q(z ,R)). (22) 

Remark 3.1 The function p R satisfies |7p with the same constants 

Henceforth the symbol C will denote a positive, finite constant that may 
vary from line to line; the relevant dependencies on the data will be specified 
whenever it will be required. Notice that the dependence never occurs on 
the unknown functions u, p or 9. 
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4 Proof of the auxiliary results 

4.1 Proof of Lemma [372] 

1) Considering the change of variables 

Q(e , r) = B(y , r) x (s - r 2 , s ) -> Q(*o, i2) = S(z , x {h - R 2 , to) 

{y, s) i-> ( x + — (y-y ),t + ( — J (s-so) 



r 



we get (fl8]) - (l20]) : observing that 

<9 s ii K + (u fl ■ V„)u B + V y p R = (j^j {d t u + (u • V)u + Vp); 

div^^A,"*) = ^ dxv y (ji{6)D v u) = (^j div(/i(0)£>u); 



» R (9 R )\D y u R \ 2 = (f \ M\Du\ 2 . 

2) Considering the change of variables 
Q(z ,R) -> Q(eo,r) 

OM) ^ H/o + -^(ar - ar ), «o + (-^) (* - *o) J 
the Jacobian is (r/R) n+2 . Then, for every m > 1 and it G L™((5t), 

/ |u(a; + -(y-y ),t + (-) (s - s ))\ m de = I \u\ m dz. 

JQ(e ,r) r \ r / JQ{zo,R) 

In particular, it follows 

r(/ \u R \ a ) 1/a = R(i \u\ a ) 1/a 

J Q(e ,r) J Q(z ,R) 

l p R\(n+2)/n\n/(n+2) _ R 2/[ |^| (n+2)/n\n/(n+2) 

Q(e ,r) JQ(z ,R) 



r 2 



r 2 



jQ(e ,r) JQ{z ,R) 
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Consequently, (|2"Tj) holds. Also we have 

Thus we obtain (1221) and Lemma 13.21 is proved. 

4.2 Proof of the decay estimate (Theorem 13.11) 

Suppose the opposite, i.e., there exist a sequence {(e m ,7 m )} and solutions 
9 m ) of ©-© in fi m x]0,T m [ and Q(z m ,R m ) CC fi m x]0,T m [ such 

that 

-R m |(u m )Q( ZmiJ R m )| < 1, -R m |(6 , m)Q(^ m , J R m )| < 7m/4 (23) 

e m -> 0, /^(zmj u m ,p m , m ) + d m R m = 7 m -> 0, (m -> oo) (24) 

Yc;R m (z m ;u m ,p m ,9 m ) > C^ a -f m , (25) 

where d m = c\(f m ] Q m x]0,T m [). Applying the change of variables 

Q = Q(0, 1) <-> Q m = Q(z m ,R m ), 
e=(y,s) 2 = (ac,t) 

we introduce the transported functions: 

w m (e) := — (u m (z) - (u m ) Qm ), 

7Tm(e) := —(Pm(z)-(p m ) Bm ), B m = B(x m ,R m ), 

7m 

x m (e) := — (9 m (z) - (#m)oJ. 

7m 

These functions satisfy the system (in the sense of distributions) in Q 

d s w m + R m ((7 m w m + a m ) • V y ) w m - div^(/i(^-x m + b m )D y w m ) = 
= -V^7r m + g m , div^w m = 0, 
d s H m + -R m (7 m w m + a m ) ■ V y x m - kA 

y^m 

= RmlmK^^m + b m )\D y W m \ 2 . 

(26) 

Here g m (e) = ^f m (z), a m := (u m ) Qm and b m := (6 l m ) Qm . 
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From (j2"4"|) we can extract subsequences, still denoted by the same symbols, 
such that 

w m -winL fl (Q), vr m - vr in L^/ n (Q), x m - x in L a: ' (Q) 

and taking into account that by definition (w m )g = {xm)q = and (7r m )s(s) = 
0, s G (—1, 0) and also the mean integrals remain zero for the weak limits, it 
follows 

1 > liminf Yi(0; w m ,7r m ,x m ) > Yi(0; w, 7r, x); (27) 

m— *oo 

liminf y,(0; w m , vr m , x m ) > C*q a . (28) 

m^oo 

In order to obtain more information about (w, 7r, x) we wish to pass to the 
limit the system satisfied by (w m , ir m , x m ) considering the weak formulation 

- / w m -d s vde+ / /i(^-x m + b m )D y w m : D y wde = 

JQ JQ Rm 

= Rm I (7m W m + 

■/ Q ^ Q JQ 

for all v G C£°(Q); (29) 

- / >c m d s (t)de-k I x m A y (f)de = R m / x m (7 m w m + a m ) ■ V^cfe + 
J Q JQ JQ 

+R mlm [ /i(^x m + b m )\D y w m \ 2 <pde } for all G C °°(Q). (30) 

JQ *Mn 

The dependence of the viscosity on the temperature does not allow to 
proceed as in the N-S system (under isothermal behaviour). In the following 
we use the local energy inequalities fl9|)- f|T0|) for (u m , p m , 6 m ) in order to obtain 
some indispensable estimates on the required solutions under some standard 
conditions on the functions involved therein. 

4.2.1 L°°(-^ 2 , 0; L 2 (5(0, ?))) n L 2 {-q 2 , 0; H^O, ?))) estimate for w m 

The local energy inequality flU]) for (u m , p m , 6 m ) can be rewriten for (w m , 7r m , x m ), 
for a.e. s g] — 1, 0[, 

[ l Wm ^' S ^ ip 2 (y,s)dy+ I I fi(j^-x m + b m )\D y w m \ 2 (p 2 dy dr < 

J B 2 J-xJb Rm 
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< 



R m / / M?"^ + b m )ipD y w m : (w m ® V y (p)dydr + 

+ y y |w m |V(3><£ + RmilmWm + 3m) * V y ip)dydr + 



+2 / / -K m (pw m -V y (fdydT + / / g m -w m (p dydr, V(p G C™(Q), (31) 

■7-1 7-1 JB 

where = 5(0, 1) C IR n . Using the Cauchy-Schwarz and Young inequalities 
we deduce 



/ l Wm ^' S )l ip 2 (y,s)dy + ~ [ [ {i(j^x m + b m )\D y w m \ 2 ip 2 dydT < 

-of I ^(^-^m + b m )\w m (g)V y (f\ 2 dydT + 

* J-l J B -Km 



fS 

+ 



J J \w m \ 2 ip(d s ip + R m (-f m w m + a m ) ■ V y ip)dydr + 
+2 / n m ipw m ■ Vyipdydr + / g m • w m ip 2 dydr. (32) 

J-l J B J-l J B 

Thus, taking a conveniently chosen tp and using (127)) under (1T2")) and also (12"31) 
we obtain 

ess sup ||w m (s)|| 2) B(o,?) + ||V y w TO || 2 ,Q(o, t ) < C+ ||g m ||2,Q, (33) 

se]-? 2 ,o[ 

with C = C{fj,#,n*). 

4.2.2 L°°(-q 2 , 0; ^(5(0, ?))) estimate for + 6 m and x m 

After letting ( — > + , the local energy inequality (ITUl) for (u m ,p m , 9 m ) can be 
rewriten for (w m , n m , x m ), for a.e. s e] — 1, 0[, 

/ (jT Xm + ^( y,S ^ dy ~ lm j j (lm^m + Rmbm)^m-V y 1pdydT + 

+ / / ( jjr K ™ + b m ) (d s i/j + kAyijj + am ■ V y V) dydr + 



'-1 J B \ 1L m 




+7m / / (J>{^>(rn + b rn )\D y w rn \ 2 'if)dydT. 



1 
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Thus, taking a conveniently chosen ip and using the Gronwall Lemma, (I33|) . 
(J2ZD, ([23D and 0, we obtain 



ess sup ||— Xm + b m \\i : B(o s ) < l m C(l + ||g m || 2 ,Q) exp[C + i2 m 7 m ||g m || 2 ,Q]. 

se]-? 2 ,o[ %i 

(34) 

By definition ffTTT) we have 



1/2 d2 / /• \ 1/2 

gm |2 de ) = _Z» II ^2 dz 

lm \J Q(z m ,R m ) 



p2 td/3 

< fhz dmR ^ = fhn dmR ^P < ijA-/9 __, o, asm^O (/3<A). (35) 
Finally, since 6 m > and using ( 1351) . from ( 1341) we find 

ess sup ||x m ||l )B ( , r ) < CR m J m . (36) 

se]-r 2 ,o[ 

4.2.3 Passage to the limit asm^oo 

In order to be in conditions to pass to the limit in fl2§|) - fl3"U|) . let us state 
the following convergences. From (1231) we get i? m a m — > a in R™. Moreover, 
|a| < 1 arises. Inserting (135!) into (1331 we get 



DyW m D y w in L 2 (<2(0, q)), 

and consequently 

Rmlm /J,(-^X m + b m )\D y W m \ 2 (/)de < Rm'Jm^WDyW 
JO *hn 



<'i|l2,QH"l x < 



< CRrrCim -> 0. 

Inserting (|35|) in (134")) and (1361) as tending m to infinity, we obtain 

^-x m + 6 m ->0, x m ->0 in ^(^(O,?)) (37) 

and consequently /i(^-x m + 6 m ) — > /i(0) in BL 

Then we can pass to the limit in (T29]) - (T30|) with x = and (w,7r) is a 
solution to the time dependent Stokes system (in the sense of distributions) 



d s w — jj,(0)A y w + V y 7T = — (a • V y )w = — V y • (a ® w) 



div y w = 
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in Q. 



The classical theory for the nonstationary Stokes equations [2H] claims 
from (a-V„)w G L 2 (Q) that w G W^Q) ^ L 2 ( n+2 )/ n (-l, 0; W 1 ' 2 ( n+2 )/ n ( J B)). 
Applying the bootstrap argument from (a • V y )w G h 2( - n+2 ^ n (Q) it follows 
that w G W^J H . 2)/ri (Q) C°' a (Q) for < a < 2 - n/2. Then w is Holder 
continuous in the closure of the cylinder Q(0, and the following estimate 
holds 

y ? (0;w,0,0) < C**<^ (38) 
with some constant C** depending only on /i# and /x*. 

Remark 4.1 TTie above bootstrap argument can still be applied if we use 
first the embedding U ■=— > L 2 ^ n+1 - )//n ((5) and next the L p -theory for the Stokes 
equation with RHS in the divergence form fl^j , that is, from w G \J 2<Jl+l ^ n (Q) 
it follows that w G L 2 ( n+1 V n (-l, o ; W 1 ' 2 *"* 1 ""^)). 

On the other hand, we need to extract subsequences which converge 
strongly in order to pass to the limit in the integral Y q (0; w m , n m , x m ). From 
d2U) and using ©, ([2D, dSSD, dSSD, USD, we derive that {<9 s w m } meN is 
bounded in L( n+2 )/ n (— 1 ; 0; W -1 '( n+2 )/ 2 (.B)) since the following estimate holds 

/ W m ■ d s \de < (/i # |p„W m || 2 ,g + ^mTm||w m ||2( n+2 )/n,Q+ 
JQ 

+i? m |a m | ||w m || 2 (n+2)/n,Q + W^mW (n+2)/n,Q) || Vj/V || ( n+2 )/2,Q + || gm || 2,Q || v || 2,Q; 

for every v G I>+ 2 )/ 2 (-l, 0; Wj ,(n+2)/2 (S)). Thanks to JSJ) and ([35]), the se- 
quence {w m } is bounded in L°°(-q 2 , 0; L 2 (£(0, ?)))nL 2 (-<; 2 , 0; H 1 (S(0, ?))) 
L 2 ( n+2 )/ n (Q(0,<;)). Using H\B) L 5 (B) and a compactness result [28] 

we obtain 

w m ^w inL a (Q((U))- (39) 
Therefore, using fl3"9l . f[3"5]) and f[3"oT) we obtain 

limsup y f (0; w m , vr m , x m ) < 

m— >oo 

n/(n+2) 



< F f (0;w,0,0) + limsup J -j Q{Q . |7T m - (7r m ) B(0if) |( n+2 )/™de 
< CU Q + 4 limsup ? ( ^ , |vr m |( n + 2 )/^e) 



< 



n/(n+2) 



In order to estimate the pressure, the idea is to decompose the pressure 
into two parts, which one part is harmonic and hence is smooth (see |15j). 
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Choosing v(y,s) = x(s)V</?(y), where x ^ Co°(0> 1) and G C™(B), as a 
test function in remarking that divw m = and 

— [ w m - Wip-^dyds = f V • w m ip^-dyds = 0, 
Jq as 7q as 

we obtain, for a.e. s G] — 1, 0[, 
/ M7T~ x m + b m )D y w m : D y Wyipdy - R m lm / w m <g> w m : Vj^dy = 

-Km JB 

= / 7r m Aytpdy + g m • Vy(pdy. 

JB J B 

Let g m be defined by 
/ /i(^x m + b m )D y w m : DyWyipdy - i? m 7m / w m <g> w m : Vj^dy = 

= / q m A y ifdy + / g m • Vyipdy, (40) 

which is valid for a.e. s g] — 1, 0[ and for all ip G Wq (B), and 

(vr m - q m , \ip) = 0. (41) 

In ( 14TT) . if we consider as a solution of the Laplace problem 

AyLp(s) = \(7T m - g m )(s)| 2/n sign((7r m - q m ){s)) in 73, 

we obtain n m = g m a.e. in Q. In ( HOI) , if we consider <£> as the unique solution 
of the Dirichlet-Laplace problem 

Ayip(s) = \q m (s)\ 2/n sign(q m (s)) in B(0,q); 

<p(s) = on B\B(0,q), 

then using ((Tj) and the Holder inequality, it results 

/ |g m | (n+2)/n < ||gm||2,B(0,,)l|V^I|2,B(0,c) + 
JB(0s) 

+ (v*\\VyW m \\ (n+2 )/n,B(0,<;) + Rmlm || W m ||2(n+2)/n,B(0,f)) II ^Vll (n+2)/2,B(0,g) ■ 
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It is known that the following estimate holds 

l|V v ¥>(s)||( n +2)/2,B(0,0 + l|V^(s)|| (n+2 )/2,B(0,c) < C \\q m (s) || (i+ 2 ) /rii B(0, ? ) ' 

for some constant C > 0. So it follows 

ll 7r m||(n+2)/n,_B(0,f) — ^ (^m A ^ + || V y W m || ( n +2)/n,B(0,f ) + Rmlm) , 



taking into account ( 1351) and (1241) . Applying the Holder inequality, we con- 
clude that 

n + 2 

limsup/ |vr m | (n+2)/n c/e < Cilimsup ( / |V„w m | 2 de) , (42) 

m-»oc jQ(0, c ) nwoo VQ(0,?) / 

with d = 

Taking in (jB2D V 9 £ Cg°(Q(0 s 2?)) the cut-off function such that 93 
I in Q(0, f ) and |V^| < C/^ 2Q /( n+2 \ |%>| < C/^ a ^ n+2 ^ inQ(0,2?), it results 

/!# / |V^w m | 2 rfyrfg < 4 ? +2) / |w m | 2 cfr/ds + 

+ , 2 x / \ir m w m \dy ds + / |g m | 2 d?/ds. 

Using (135]) . the Holder inequality and observing that 

(n+2 
— -J^vT) / k«w»|de ) " < ^ / |vr m | (n+2)/n rfe + 
m# ^ n+2 > y Q(0 ,2 ? ) / 2d j Q{0>2q) 



c 



r 2a/n \ 1 
S 1 VQ(0,2 S ) 



|w m | a rfe 



n + 2 

an 



n+2 _ v n+2 

a. 



2j„ / ° 2 



we conclude 

limsup^y |V y w m | z <ieJ < lim sup ( -f |w m | a <ie 

2C*i m-»oo JQ(0,2<;) 



Q(0,c) / s m ^°° VQ(0,2?) 

1 /" n + 2 

H — — lim sup -f |7r m | » ae, (43) 



with C 2 = C(/z#,/i#). 
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Introducing (T431 into (jl2l and applying (1551) . we obtain 
lim sup 



/ |7r m |( n+2 )/^e < tf 3 (2c) a + ilimsup/ |7r m |( n+2 ^e, 



with c 3 = c 1 c 2 2 2Q /«d: +2)/n . 

Iterating over i = 1, • • • , k, where k is such that 1/2 < 2 k <; < 1 we derive 

k 

lim sup-/ \ix m \^de < 4C 3 ^2 i(a ~ 2 V + lim sup -f |7r m |^rfe. 

m-»oo JQ(0s) i = i m->-oo jQ(0,2 fe <;) 

Considering that 

/ \-K m \ (n+2)/n de<C2 n+2 ( |7r m | (n+2)/n rfe < C2 n+2 

jQ(0,2 k <;) JQ 

it results 

hmsup f K^^de < 3 / + C < Ctt 2)/n - 

Hence we obtain 

lim sup F c (0; w m , 7r m , x m ) < C«^ Q + C,** 1 ? < (C** + C***)^. 

m— »oo 

This is a contradiction with (|28|) if we set C* > C** + C***. Theorem 13.11 is 
proved. 

4.3 Proof of Lemma 13.11 

(i) We prove (fT6l) by induction on m. The proof of f|T6|) when m — 1 follows 
from Theorem I3.1[ observing that (1T41) holds and 

CU a = cu a ~V < s s . 

Now we suppose that ( |T4l) holds and, for z = 1, • • • , m, we have 

Y^ R (z;u,p,9) < q iS (Y R (z;u,p,6) + c x (f)RP). (44) 
Thus (glD and flBD imply that 

F^(z;u,p^) + c A (f)(^ < ^(F i? (^u,p^)+c A (f) J R' 3 ) + c A (f)^ 

< g iS Y R (z;n,p,9) + c x (f)R^(<; iS + ^) < 7 , 
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considering that (TT5"|) holds. 

We can apply Theorem 13.11 and subsequently (jUJ) resulting 

Y, i+iR (z;u,p,9) = Y^ iR) (z;u,p,9)<C^(YM^,P,0) + cx(i)(^) 

< C^ a (<; i5 (Y R (z; u,p, 9) + c A (f)i^) + c A (f) 

< ^ 5 ( ? «Y fl (z;u,p,0) +2^c A (f)^) 

< <; {i+1)5 (Y R (z;u,p,9)+c x (f)RP), 

which completes the proof of (i) in Lemma 13.11 

(ii) Let p e]0,sR] be arbitrary and mGNbe such that 

m+l ^ P_ ^ m 

Proceeding as in [TJ] we have 

Y p (z;u,p,0) < C(<;)Y qmR (z;u,p,9) < (see 



< C7(?)(— ) (y*(z;u,p,0) + c A (f)i^ 



< ^)(j£) (^;u,p,0) + c A (f)iO, 
which concludes the proof of Lemma 13.11 



5 Proof of Theorem 12.2 



Choose c such that (fT5l) is fulfilled. For instance, taking 5 = a/2 and (3 — A/2 
it follows 

? a/2 < ^ ? a/2 + ^A/2 < L 

Thanks to Theorem 13.11 with the above chosen ^ and /?, there exist e = 
£((7*, a, A) and 7 = 7(C*, a, A) such that for any proper weak solution (u, p, 9) 
of the N-S-F problem in Qt satisfying 

Q(z,R)ccQ T , 0<R<e, 
R\(u) Q{2;R) \<l, R\(9) Q(Z)R) \<j/4 } (45) 
Y R (z; u,p,9) + c x (i)R x / 2 <7 
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and Lemma 13.11 can be applied. Indeed, for all < p < $R we get 

Y p (z; u,p, 6)<C (£f' 2 (Y R (z; u,p, 9) + c A (f)i? A / 2 ). (46) 

Let A = £7/8 be the desired constant. By hypothesis ffl3|) for each 
z o £ Qt there exists a constant Rq > that can be chosen such that Rq < e/2 
and for all R < Rq we have 



g(2b,i2) cc g r; it;(y^ ;u, P ,e) + c A (f)i? A / 2 ) < -J- + -J- = -±. 

By the continuity of z 1— > ^(z; u,p, 6 1 ) at zo> there exists a neighbourhood 
of Zo, O(zq), such that for all z G 0(z o ) we have 

R(Y R (z;u,p,9) + c x (f)R x / 2 ) < e 7 /4. (47) 

Indeed, there exists < i?i < _R such that for all R < R\ and Q(z, R) C 
0(zo) the relation (1471) is satisfied. Applying Lemma T3.2I with e = (y,s) and 
r = e in (l2"Tj) - (l2"2"j) and using (T4T1) we obtain 

Y r (e;u R ,p R ,6 R ) + c x (f R ;Q r )r x / 2 < 
<jh R (z;u,p,e)+^^I^c x {t)j < 7 /4 (J2<r), 

with Q r denoting the transported domain. Let us consider the transported 
solution (u R ,p R ,6 R ) and < r < min{e, 4/7}. Then the assumption (T4"5]) is 
fulfilled regarding that r < e, and 

r|(u R ) Q(e)r) | < ry r (e;u* p*,0 fl ) < 1, 
r|(^ (e>r) |<r r (e;u ii ,p«0<7/4, 
y r (e; u fi ,p R , ^) + c A (f*; Q r )r A / 2 < 4? r (e; u* 0*) + c x (i R ; Q r )r x ' 2 < 7. 

Now using (146]) it results 

Thus u is Holder continuous with exponent a/4 in a neighbourhood of 
e, taking into account the parabolic version of the Campanato criterion [4j 
Theorem 1.2]. Theorem 12.21 is proved. 
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6 Proof of Theorem 12.3 



We begin by recalling the following result which plays a central role in the 
proof. 

Lemma 6.1 f[T2l Lemma 11]) Let f E L\ oc {Q T ) and, for0<d<n + 2, 
let 

F = {z E Q T : lim sup Br* i \f\dxdt>0}. 

R-> 0+ JQ{z,R) 

Then, we have H d (F) = 0. 

Let z E S be arbitrary and, for B > 0, denote 

A = Ar n [ l u l a > B = — W [ \p\^, D = Ar n l H a '> 

^n-K JQ(z,R) ^n-tl J Q(z,R) ^n-K JQ(z,R) 

with u n denoting the measure of the unit n-dimensional ball. 
First, we observe that 

5 

Sc\Jf u 

where each is the set of points of S that verify the following i-case. 

1. [A < 1, B < 1, £) < 1] This case is impossible due to 

< A < lim sup B{B~ 2/a + B(BT 2n/{n+2) + BT 2/a ')) 

R^0+ 

= lim sup(i2 1-2/a + B 4/{n+2) + B 2/a ) = 0, 

R^0+ 

taking a > 2 into account (cf. f[T2"j) ). Then we find that Fi = 0. 

2. [A < 1, 5 < 1, D > 1] Since A a '/ a < 1 and _B a V( n + 2 ) < l, W e get 

< A a ' < 4 lim sup (B (1 - 2/a)a ' + B ia ' /{n+2) + B 2a ' ,a D) 

R 2a'/a r 

< 41imsup — / \9\ a dxdt. 

U n B n jQ{z,R) 

Then from Lemma O it follows n n ~ 2/{ - a ~ 1) {F 2 ) = 0. 
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3. [A < 1, B > 1, WD] Considering that 

2(n + 2) , 3 n + 2 2(n + 2) 

1 < -i - 1 < a < - < < 3 < a < -i '- 

n + 4 2 n n 

implies that B«V("+2) < B, we get 

< A"' < 4 lim SUp (R^- 2 /^' + R^'/(n+2) B + ^a'/a^, 



R 2a' /a 



n + 2 



< 41imsup / |p|~ + \Q\ a dxdt. 

R^0+ U n R n Jq( Z) R) 

Then from Lemma O it follows / H n ~ 2 ^ a ~ 1 \F 3 ) = 0. 

4. [A > 1, 5 < 1, VD] Arguing as in the above two cases, now considering 
that A a '/ a < A, we get 

< A a ' < 4 lim sup (R (1 - 2/a)a 'A + # 4a '/(«+2) + R ^'h D 

R (l-2/a)a' r 

< 41imsup — / |u| a + \9\ a dxdt. 

h-»o+ u n R n Jq(z,r) 

Then from Lemma O it follows ?f-(a-2)/(a-i) (i? 4 ) = o. 

5. L4 > 1, 5 > 1, VD] Analogously to the above cases, now with A a '/ a < 
A and 5°V(«+2) < B, we get 



< A a ' < 4 lim sup f i? (1 - 2/a)a ' A + R^'/^+^B + i? 2 ^/)' 

R (l-2/a)a> r . n+2 , 

< 41imsup — — / (|u| a + + \6\ a ) cfacft. 

_R^0+ U n R n jQ(z,R) ^ ' 



'Q(z,R) 

Then from Lemma O it follows ft»-(«-2)/(o-i) (_p 6 ) = 0. 
Next, from the properties of measure, we have 

5 5 
i=2 i=2 

supposing that d = maxjrfj, i = 2, • ■ ■ , 5}. 

Then, by the assumption (fl~2]) we have a < 4 (n = 2, 3) and we conclude 
that 7-t n ~( a ~ 2 y( a ~ 1 ) (S) = 0. Therefore by definition of parabolic Hausdorff 
dimension the proof of Theorem 12.31 is finished. 
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